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The unprecedented light curves of the Kepler space telescope document how the brightness of 
some stars pulsates at primary and secondary frequencies whose ratios are near the golden mean, the 
most irrational number. A nonlinear dynamical system driven by an irrational ratio of frequencies 
generically exhibits a strange but nonchaotic attractor. For Kepler's “golden” stars, we present 
evidence of the first observation of strange nonchaotic dynamics in nature outside the laboratory. 

This discovery could aid the classification and detailed modeling of variable stars. 

PACS numbers: 05.45.Tp, 05.45.Df, 97.10.Sj, 95.75.Wx 


As a byproduct of its vastly successful search for ex¬ 
oplanets [lJi3, the Kepler spacecraft has revolutionized 
stellar photometry by its high precision tracking of the 
brightness of 150 000 stars nearly continuously for four 
years. Both the quantity and quality of the Kepler data 
enable new stellar discoveries. The Kepler stars include 
cosmic distance markers (or standard candles) such as a 
Cepheid and 41 RR Lyrae variable stars. While some of 
these stars pulsate with a single frequency, Kepler ob¬ 
servations confirm that others pulsate with multiple fre¬ 
quencies. Several of these stars, including the RRc Lyrae 
star KIC 5520878, pulsate with two principal frequencies, 
which are nearly in the golden ratio [o]. As the most ir¬ 
rational number, the golden ratio can have significant 
dynamical consequences. For example, by the KAM the¬ 
orem, dynamics with two frequencies in the golden ratio 
maximally resist perturbations @]. Furthermore, nonlin¬ 
ear systems driven by two incommensurate frequencies, 
forming an irrational ratio Q , and especially the golden 
ratio, exhibit distinctive dynamics that exist between 
order and chaos. Strange (because fractal) nonchaotic 
attractors characterize this distinctive dynamics 00. 
Here we report evidence for strange nonchaotic behav¬ 
ior in the pulsations of stars like KIC 5520878, making 
this golden star the prototype of a new class. Strange 
nonchaotic attractors have only been seen in laboratory 
experiments but never before in nature. The observa¬ 
tion of stellar strange nonchaotic dynamics provides a 
new window into variable stars that can improve their 
classification and refine the physical modeling of their 
interiors Q 

Periodically driven damped nonlinear systems can ex¬ 
hibit complex dynamics characterized by strange chaotic 
attractors, where strange refers to the fractal geome¬ 
try of the attractor and chaotic refers to the exponen¬ 
tial sensitivity of orbits on the attractor. Quasiperiodi- 
cally driven systems forced by incommensurate frequen¬ 
cies are natural extensions of periodically driven ones 
and are phenomenologically richer. In addition to pe¬ 
riodic or quasiperiodic motion, they can exhibit chaotic 


or nonchaotic motion on strange attractors. Although 
quasiperiodic forcing is not necessary for strange non¬ 
chaotic dynamics, the first experiment to demonstrate 
a strange nonchaotic attractor [12] involved the buck- 
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FIG. 1. Spectral content, (a) Raster pattern of disks with 
brightness proportional to the KIC 5520878 stellar flux at 
intervals of St ~ 1766 s ss 0.5 h. (b) Fourier transform mag¬ 
nitude of stellar flux sampled at intervals St has primary and 
secondary frequencies (red arrows) at /i ~ 1/(0.266 d) and 
/2 ss 1/(0.169 d), where / 2//1 ~ 1.58 r; ip is the golden ratio. 
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FIG. 2. Attractor reconstruction. Two views of a three-dimensional plot of normalized flux Fn at t = nSt successively delayed 
by r = 2<5t. This delay coordinate embedding “unfolds” the time series into a warped torus, suggesting two-frequency nonlinear 
dynamics. Equal-sized spheres locate data. Rainbow colors code time, from red to violet. Flux triplets straddling data gaps 
appear far from torus. 


ling of a magnetoelastic ribbon driven quasiperiodically 
by two incommensurate frequencies in the golden ratio 
/ 1//2 = ip = (1 + a/ 5)/2 ~ 1.62, which is the irra¬ 
tional number least well approximated by rational num¬ 
bers, as its continued fraction expansion tp = 1 + 1/tp = 
1 + 1/(1 + 1/(H-)) is all Is. Here we apply the time se¬ 

ries analysis techniques of the ribbon experiment to ana¬ 
lyze the dynamics of several Kepler multi-frequency vari¬ 
able stars, including the RRc Lyrae star KIC 5520878, a 
blue-white star 16 000 light years from Earth in the con¬ 
stellation Lyra whose brightness varies as in Fig. [IJa). 

Most classical Cepheids and RR Lyrae are single-mode 
periodic and radial pulsators. For many years, the only 
known exception was long-term Blazhko light curve mod¬ 
ulation found among some RR Lyrae [i]| . Recently, bet¬ 
ter quality photometry has uncovered various multi-mode 
pulsations within these variable stars, possibly from non- 
radial modes. A few dozen such stars have been docu¬ 
mented 0 for mysterious behavior involving frequency 
ratios near 0.62, apparently without connection to the 
golden ratio p ~ 1/0.62 or the possibility of strange non- 
chaotic behavior. We refer to the former stars as golden 
stars and to the latter stars as strange nanchaotic stars, 
and we recognize that these classes may overlap. 

We downloaded the Kepler Input Catalogue 5520878 
data from the Mikulski Archive for Space Telescopes. 
The light curve, or flux F (in detector electrons per sec¬ 
ond) versus time t (in days), was sampled every 5t = 1766 
seconds or about once each half hour, for most of four 
years, including over 5000 primary periods. To keep its 
solar cells in sunlight, the Kepler spacecraft rotated 90° 
every quarter solar orbit causing the star to illuminate 
different pixels of its charge-coupled device. Different 
pixels have different responses and appear to have con¬ 
tributed to shifts and skews in the data segments. In 


addition, events such as cosmic ray strikes, monthly sci¬ 
ence data downlinks, and safe modes caused both small 
and large gaps in the data. The best detrending practice 
for planet searches jl5| cannot be employed for variable 
stars, such as KIC 5520878, because each pulsation is in¬ 
trinsically different. Instead, we manually processed the 
raw time series F by detrending and rescaling the data 
segments to zero mean and unit variance to obtain a re¬ 
fined or normalized time series iby. 


A natural strategy to test the flux time series for 
strange nonchaotic behavior is to estimate both the di¬ 
mension of the underlying attractor and the largest Lya¬ 
punov exponent of the orbits. A fractional dimension 
would imply fractal geometric scaling and strangeness, 
while a positive Lyapunov exponent would imply the 
divergence of nearby orbits and chaos 7]. Although 
these metrics are difficult to estimate [16J, we investi¬ 
gated both. 


For the fractal geometry 0 , we calculated the infor¬ 
mation dimension, which typically delivers the most ac¬ 
curate dimension estimates for short noisy time series 181 
by quantifying how the information needed to specify a 
set of points scales with the area that contains them. To 
estimate it, we partitioned a section of the attractor into 
many tiny boxes, computed the fraction of points in each 
box, and calculated the corresponding information or en¬ 
tropy. The logarithmic scaling of information with slope 
d ~ 1.8 in the range 1 < d < 2 is consistent with some 
fractal structure. (In contrast, for pure quasiperiodic mo¬ 
tion, the section is a sinusoidal curve of dimension d = 1, 
and for completely noisy data, the section is a uniform 
area of dimension d = 2.) 


For the largest Lyapunov exponent 11711. we imple¬ 
mented both the Wolf [19| and Rosenstein [20] algorithms 
and found the latter more accurate and stable when 
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FIG. 3. Spectral scaling, (a) Fourier transform magnitude of stellar flux strobed at the primary period Ti has a wide distribution 
of peaks. Peaks above threshold a (horizontal line) are highlighted (gold), (b) Log-log plot of number of super-threshold peaks 
N versus the threshold a has a linear regime (in leg below knee of curve) that indicates power law scaling and fractal structure. 
Spectral exponent slope —2 < s < — 1 indicates strange nonchaotic dynamics. Deviation for low thresholds arises from light 
curve’s finite temporal resolution. 


tested on known time series. Beginning with the Fig. [2] 
attractor reconstruction, we found the nearest neighbors 
of all points not too close together. We then computed 
the average separation of all neighboring points as a func¬ 
tion of iteration. A linear increase in the logarithm of the 
average separation versus iteration number would have 
indicated a positive largest Lyapunov exponent. How¬ 
ever, we did not observe such a linearity, which is consis¬ 
tent with a zero largest Lyapunov exponent [20I ] . 

Since these two obvious metrics are difficult to esti¬ 
mate, we focused instead on the spectral scaling of the 
dynamics, which is a more reliable measure that has tra¬ 
ditionally been used to identify strange nonchaotic mo¬ 
tion 13, 21, 22| ■ We analyzed the spectral scaling in two 
independent data pipelines with comparable results. One 
pipeline used Mathematicci 10 software while the other 
used MATLAB, Period04, and custom C++ software. 

From a Fourier transform Fn [/] of the normalized time 
series F,y [t] • we identified the two most significant fre¬ 
quencies in the light curve, as shown in Fig. [lib), where 
red arrows indicate primary and secondary frequencies 
at /1 ~ 1/(0.266 d) and f 2 ~ 1/(0.169 d), correspond¬ 
ing to primary and secondary periods of T\ ss 6.41 h 
and T 2 ~ 4.05 h, where / 1//2 = T 2 /T\ ss 1.58 ~ ip is 
nearly the incommensurate golden ratio, in agreement 
with previous work [3l23[. (Other spectral peaks corre¬ 
spond to harmonics of f± and linear combinations of /1 
and f 2 .) We next plotted three-dimensional delay coor¬ 
dinates {+/v[t], F)v[t — r], Fjv[t—2r]} at times t = ndt and 
delay r = 2 St to reconstruct the time series’ underlying 
attractor Q. The resulting Fig. [2] torus is consistent 
with quasiperiodic forcing of the dynamics. Points far 
from the torus straddle gaps in the data, and the warp of 
the torus reflects the nonlinearity of the light curve [l7(. 

One way to take a Poincare section of a state space flow 
is to intercept it with a plane, but a second way is to sam¬ 
ple or strobe it at a fixed frequency. As in the original 


strange nonchaotic experiment 12], we strobed the data 
at the primary period T\ and wrapped it by the secondary 
period T 2 to produce the section. To qualitatively check 
the section dynamics 0 , we tracked the normalized flux 
differences A Fn = FN[(ni+n)Ti] — FN[(n 2 +n)Ti] of two 
points m and n 2 as a function of iteration n. For many 
steps, this difference is small, briefly increases to large, 
and rapidly returns to small, a characteristic nonchaotic 
behavior. (In contrast, for strange chaotic attractors, the 
difference would start small, increase to large, and there¬ 
after remain uncorrelated, without the repeated rapid 
close returns.) 

To quantify the section geometry, we examined the 
magnitude of the discrete Fourier transform of the sec¬ 
tion (or strobed) data, as in Fig. [Ha), and recorded the 
number of peaks N above a threshold er. The power law 
scaling of this spectral distribution N = Nocr s , over 2 
decades of thresholds and 3 decades of number, with an 
exponent s ~ —1.5 in the range —2 < s < —1, as in 
the Fig.[3lb) log-log plot, is a classic signature of strange 
nonchaotic dynamics [9| [l0|. The power law means that 
the distribution of peaks is scale free, with a range of 
large and small peaks, presumably reflecting the fractal 
nature of the underlying strange attractor [2J, |25( ■ The 
deviation at low thresholds results from the light curve’s 
finite frequency resolution. 

The cumulative case for strange nonchaotic dynamics 
includes the strobed or section spectrum, the spectral 
exponent, the distinctive orbit separation behavior, and 
to a lesser extent the section dimension and Lyapunov 
exponent. The identification is challenging because of 
inevitable noise, limited and imperfect (despite unprece¬ 
dented) data, and the smallness of the secondary fre¬ 
quency amplitude compared to the primary frequency 
amplitude (as suggested by the Fig. [Hb) spectral peak 
heights). Quasi-periodic dynamics naturally flows on a 
two-dimensional torus, with the two frequencies driving 
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FIG. 4. Other stars, (a) Spectral scalings of the four Ke¬ 
pler golden stars all exhibit power law scaling, over 2 decades 
of thresholds, and nearly collapse onto each other when the 
axes are proportionally scaled (as denoted by primes), (b) 
By contrast, spectral scalings of two Kepler non-golden stars 
with primary and secondary frequency ratios near 3/2 do not 
exhibit power law scaling but are in fact concave or multi¬ 
scaled. 


the toroidal and poloidal motions; quasi-periodic forcing 
of a nonlinear system naturally inhabits a wrinkled torus 
with fractal cross section, but in this case the wrinkles 
may be small (or the strangeness higher dimensional), 
and standard analysis might miss them. 

We tested our time series analysis by analyzing a va¬ 
riety of null hypotheses that generated surrogate data 
sets of artificial light curves of flux versus time, including 
both nonparametric and parametric models. Nonpara- 
metric models included a classic phase randomization of 
the original time series, which was the inverse Fourier 
transform of the phase randomized Fourier transform of 
the light curve [26|. The strobed section and spectral 
exponent easily discriminated between the surrogate and 
original data. 

Parametric models included an ideal quasiperiodic 


curve having either the KIC 5520878 near-golden fre¬ 
quency ratio or the ideal golden frequency ratio, with 
or without noise added to the fluxes or the times. The 
spectral distributions of these artificial data sets were not 
strange nonchaotic, demonstrating that the KIC 5520878 
pulsations are not simply noisy quasiperiodic. To check 
the effects of the many small and the few large data gaps, 
we linearly or sinusoidally interpolated to close them. 
This produced some artifacts in the section plot, but 
still produced a spectral exponent in the strange non¬ 
chaotic regime. Finally, we introduced the exact experi¬ 
mental gaps into the ideal quasiperiodic curve. For lim¬ 
ited ranges of additive noise and thresholds, this did pro¬ 
duce spectral power law scaling, but the strobed spectra 
remained qualitatively distinct, and strange nonchaotic 
dynamics remained the best explanation. 


We expanded this analysis to five other multi¬ 
frequency variable stars in Kepler's field of view. The 
three additional golden RRc Lyrae stars KIC 4064484, 
KIC 8832417 and KIC 9453114, whose frequency ratios 
well approximate the golden mean, all exhibit signatures 
of strange nonchaotic dynamics, and their spectral scal¬ 
ings nearly collapse onto each other when the axes are 
proportionally scaled, as in Fig. Oja). In contrast, two 
non-golden RRab Lyrae stars KIC 4484128 and KIC 
7505345 j27| . whose frequency ratios well approximate 
the simple fraction 3/2, exhibit qualitatively different dy¬ 
namics, as in Fig. |4jb). The possibility of other golden 
stars with strange nonchaotic behavior, discovered by on¬ 
going large sky surveys like OGLE 28j or AS AS (2j|, 
promises to further classify the RR Lyrae and Cepheid 
variable stars. 


Our model-independent nonlinear analysis of the light 
curves of such variable stars is complementary to de¬ 
tailed nonlinear hydrodynamic models of the stars them¬ 
selves [30]. The strange nonchaotic signatures of variable 
stars like KIC 5520878 may elucidate phenomena like 
the Eddington valve mechanism 3l| thought to under¬ 
lie their pulsations: variations in the opacity of the star 
might quasiperiodically modulate the normal hydrostatic 
balance between pressure outward and gravity inward, 
thereby generating light curves somewhere between order 
and chaos that the best models will need to reproduce. 


Strange nonchaotic attractors have been observed in 
laboratory experiments involving magnetoelastic rib¬ 
bons [l2|, electrochemical cells |32|], electronic cir¬ 
cuits [33|, and a neon glow discharge |34j, but never be¬ 
fore in non-experiments in nature. The pulsating star 
KIC 5520878 may be the first strange nonchaotic dy¬ 
namical system observed in the wild. 
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